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We calculate the shot noise power in a junction of a network of Majorana bound states (MBS) with 
a normal metal. These Majorana bound states are on the border of alternative ferromagnetic and 
superconducting regions at the quantum spin Hall insulator edge. We analyze different realization 
of MBS networks including a few isolated ones or in a chain allowing for the limit of weak and 
strong tunneling. The conductance, as well as the shot noise are considerably stronger than those 
of a weakly coupled normal-superconduting junction, being a hallmark of the MBS. We find that 
the Fano factor is quantized F — 2 when one MBS member of a pair is coupled to the lead, however 
if both MBS members are coupled F is non-integer. 

PACS numbers: 73.50.Td, 73.23. Ad, 73.63.-b 



Introduction: Recently discovered topological materi- 
als in combination with a superconductor can host Majo- 
rana fermions [l|-l8|]. Interest in such quasiparticles is due 
to the possibility of non-Abelian statistics which they sat- 
isfy 0. The most accessible case are Majorana Bound 
States (MBS) at the 1-dimensional edge of a quantum 
spin Hall insulator [ljj Two Majorana fermions can 
form a usual Dirac fermion. However, to detect a state of 
Majorana fermion which has no definite c harge, r equires 



non-local measurements. It was suggested [12|, ll3( that a 
tunneling probe can detect the presence of MBS. While 
the tunneling probe is a local measurement it can detect 
interference effects between various MBS. Recently an ar- 
ray of alternating ferromagnetic and superconducting re- 
gions at the QSH edge were considered 15 , 1(| . This ar- 
ray may appear due to local density fluctuations at a ran- 
dom FM/SC boundaries that host MBS. The tunneling 
characteristics (conductance and tunneling current) were 
calculated for a network of coupled Majorana fermions. 
Both isolated [l6| and random chains [l5[ cases were an- 
alyzed. Here we show how Majorana fermions can be 
detected by shot noise measurements. For this we calcu- 
late the noise power for different networks of MBS and 
present a general formulation of this problem based on 
Keldysh technique. 

The Hamiltonian: We consider the tunneling between 
the normal metal lead and a realization of Id Majorana 
chain [III . EMU. The coupled Majorana state network 
is described by the Hamiltonian 



H 



M 



= 2 tijlilj 



(1) 



where 7^ are operators for the i-th Majorana fermion, 



satisfying 7$ = 



h ' 



1, and tn are elements of 



an antisymmetric matrix t for the coupling between the 
i,j MBS. Disorder may be introduced as random near- 
est neighbor coupling. The 7$ can be written as a Bo- 
golyubov transformation for the quasiparticles in the su- 
perconductor composed of creation ^(a;) and annihila- 



tion ^ a {x) operators of electrons at position x and spin 
a. The coefficients f<r,i(x) in this transformation are the 
eigenvalues of the Bogoloyubov de Gennes Hamiltonian 
with zero energy. The general form is 



dx(U ti (x)V a (x) + f:Ax)*Ux)) (2) 



The tunneling between the normal metal elec- 
trode and the superconductor is given by Ht = 
J2ka I ^ x [*fe( a; ) c lo- v ^o-( a; ) + h.c.]. Below we consider the 
energy gap A of superconductor as the biggest energy 
scale in the problem. Then for small applied voltage 
eV < A only zero-energy Majorana operator projection 
of total quasip article operator in superconductor is im- 
portant 15Ml7l| . Thus the tunneling Hamiltonian and the 
current operator in Nambu space become 



I = ^ J2 s (0)Vi7i + h.c = ie J 



(3) 



(4) 



where c(0) = (c-j-, cj,, cj, — c|) T corresponds to the normal 
lead electron operator at x = The coupling between 
the normal lead and Majorana states is given by a vec- 
tor in Nambu space Vi = (ff,i,h,i> ft,i>-f%,i) T where 
f a _i = J dxf a ^{x)t k {x) describe the interaction of MBS 
with the lead electrons and we assumed a weak momen- 
tum dependence of matrix element for tunneling to 
the lead. The Pauli r-matrices act on (c CT ; f a ,i) and (cy, 
/f f ) blocks. Total hamiltonian in addition to parts (|H3p 
includes the lead contribution which we take at voltage 
bias V. We write below the action for lead in the rotated 
Keldysh basis 



(5) 



Siead = \ dtH k c k g k c k 



where g k 1 is the inverse Green function for electrons 
in the lead. For further use we introduce the Green's 



2 



function (GF) of the lead integrated on momentum g 
■^Yj^gk that has the following Keldysh form 



q R q K 
1 g A 



(6) 



Here all entries are 4x4 diagonal matrix in Nambu 
space. In energy representation each of them has 



TfiV(O) diag(l, 1,1,1) and ig 



K _ 



a form: 

N(0)diag[tanhSi=^(l, 1,0,0) + tanh ^§^(0, 0, 1, 1)], 
where N(0) is the electron density of states in the normal 
metal lead. To calculate the noise we need to find the 
effective action as function of a quantum source. We 
combined the current operator (J) multiplied by source 
fields A=diag(Ai, — A2) on standard 1, 2 Keldysh contour. 
We obtain a combined source-tunneling contribution to 
the action 

i 

+7»V?(r,ff,-/A)c(0)] (7) 

Here the Pauli matrix a z relates quantum operators (Ma- 
jorana ji = ("fiij^fa) and lead fermions) on the two time 
contours (1,2). I is the unit matrix in the Keldysh or 
Nambu spaces. Integrating out the lead electron opera- 
tors and performing a rotation to the Keldysh form sim- 
ilar to (J6j> we arrive to a simple representation of S s . 
Keeping only the quantum source X q = (Ai — X 2 )/2 we 
obtain 



H\ q )g(T z i + ia x \ q )V 1 -i' J 



(8) 



Where 7- = (-ff,-f^) T = (jii+-f2i,-fu~j2i) T '/2. We also 
notice a subtlety: due to the sclf-conjugacy condition the 
Keldysh GF for Majorana fermions differs from ((6]) and 
is similar to that of boson fields 



G = 



G K G R 
G A 



(9) 



To zero order in tunneling to the lead the action follows 
from Eq.flTJ 



S 



dtdt'l^fGo^tt')^') (10) 



here the inverse matrix GF is 

1 



[Go 



fli-i 



id t - 2it\ 



(11) 



The total GF follows immediately if we add the zero 
source part contribution of S s to the Gq : 



G- 1 = Gq 1 - 2irV^a x gV 



(12) 



which can be used to verify the form (|9]). The effective 
action with the source term becomes 



'e// 



\ J dtj (G- 1 + Q(X C 



))7 



(13) 



Q(\) = V^b~ z (\g - Vx9°x\) + \g°x\]v (14) 

The partition function for calculations of transport and 
noise is a function of quantum field X q : Z(\ q ) = 
j D^exp[iSeff] and lnZ(A g ) = Tr ln[l + GQ{\ q )\. Tak- 
ing the variation of In Z on X q —5-0 we get expressions 
for the average current and the noise power. 



Sn(W) 



e 6 In Z e S 
2h SX q = 2hJ\~ q 



dujTr[GQ] 



2 5X q (t)8X q {t') 



hxZ 



(15) 



Before presenting explicit formulas we discuss the ap- 
proximations which are used: (a) We consider the tun- 
neling element fi a as spin independent, (b) When (|T2|) 
is inverted we neglect off diagonal fiafja 1 terms that in- 
volve oscillation with a Fermi wavevector which vanish 
upon averaging on large spacings between MBS [lil ] . (c) 
We take a wide band limit ignoring the momentum de- 
pendence of Vi. Thus we have Fy = A-KN(Q)ff5ij. (d) 
We also consider A > Tij. We proceed now to take ex- 
plicitly the trace in Keldysh and Nambu spaces leading 
to general expressions for the current and the total zero 
frequency noise power S n = Si + S2 



I 

Si 

s 2 



e 

2h 

2e^ 

h 

,2 



J duTr(TImG R ){t&nh ^ - tanh 



2T J 



(16) 



duTr(rimG R )[l - i(tanh 2 ^ + tanh 2 ^] 



— / diuTr (TReG R TReG R ) (tanh— -tanh— ) 
h J IT 2T' 



Here u± = lu ± eV, ImG R = (G R - G A )/(2i) and 
ReG R = {G R + G A )/2. The term Si mainly contributes 
to the thermal part of the noise power. In the limit of 
zero temperature (T — > 0) Si — > 0. The 52 term for 
V > T — >• defines the shot noise which is the total noise 
in this temperature limit 



8, 



hot 



reV 

/ dujTr(TReG R rReG R ) 
Jo 



(17) 



Isolated Majorana states: For single Majorana state 
the matrix i=0 and G R = G R ± = 2j{u + 2iT). First 
let us consider the zero bias conductance in the limit 
T > eV — >• 0. In this case 5^ corresponds to thermal 
noise. Directly by calculating the current (|16j) or with 



the help of fluctuation dissipation theorem: 
obtain the conductance 



2e 2 



4r 2 



h r 2 + 4r 2 



a we 



(18) 
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At zero temperature and V ^ the differential conduc- 
tance coincides with Eq. (|T8|) if we replace T by eV , 
while the shot noise for single Majorana state acquires 
the form 



le 2 T . eV 

( arctan — 

h K 2T 



2eVT 



{eV) 2 +4T 2 ' 



(19) 



More complicated formulas for S n follow when two MBS 
couple by a tunneling element t and couplings Tn, 1^2 to 
the lead. In comparison, the well known Andreev contri- 
bution to the noise of normal-superconductor (NS) junc- 
tion 18j is of order 



S 



4e 3 F r a 
h ( A> 



(20) 

or the other MBS networks 



This is much weaker then (Q2 
or by the factor (T/A) 2 <C 1. 

Fig.l presents the conductance and shot noise power as 
function U — eV/T for different realizations. With two 
coupled Majorana's the shot noise shows clear steps in 
its voltage dependence. They corresponds to the peaks in 
conductance at eV = ±2t. Therefore the measurement of 
the noise can detect Majorana fermions. In case that one 
of the links is weak: tu„k < t the noise loses its structure 
( lower right panel in Fig.l) and becomes similar to the 
case of a single Majorana (upper left panel in Fig. 1). 
In this case the sharp drop of conductance at voltage 
eV — > can be smeared by a small temperature. 

The Fano factor is defined as the ratio F = S s hot/2eI. 
For the NS junction F(V — > 0) = 2 reflecting the trans- 
mission of cooper pairs through the superconductor. The 
Fano factor corresponding to the upper right panel of Fig. 
1 (2 MBS) is plotted in Fig. 2. The zero bias value of 
Fano factor F(V = 0) = 2 does not depend on the on the 
value of t ^ and is the same integer for four MBS, the 
case in the lower left panel of Fig. 1. However, when the 
second Majorana state interacts with the normal lead, 
i.e. 7^ 0, F(V = 0) is non-integer and decreasing, as 
shown by the solid line in Fig. 2. In contrast, with odd 
number of MBS F(V -+ 0) = 0. 

Considering next the case of an infinite homogeneous 
chain, where all near neighbor couplings are identical 
U i + i = t, the Green's functions can be easily found 
[la] . The calculations of the shot noise are performed 
for t/T = 0.1 and compared with conductance in Fig. 3. 

Disordered case - isolated Majorana pairs: We consider 
next the case of a semi infinite MBS chain with a ran- 
dom ti t i+% and only the 1st MBS is coupled to the lead 
r« ~ Si t i. This system can be mapped to a system of 
independent MBS pairs [To| b y using analysis methods of 
the quantum Ising spin chain 1!| 20| • The key ingredient 
is to isolate the strongest bond at each decimation. The 
MBS pairs have an an effective hopping e™ ax ps e _e 1 
(having a typical weight) leading to peaks of conductance 
at eV — e n with weights T 2 L « e _sr ° e " /2 ; r is a start- 
ing logarithmic flow parameter chosen here as Tq = 0.1. 
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FIG. 1: (a) Upper pannels for conductance (dashed or dot- 
dashed line) and shot noise (full lines): Upper left for a single 
MB, upper right for two MBS with t = Y — upper lines (dot 
dashed and thick lines) for Tn = 21^2 = T while lower lines 
(dashed and thin lines) for Fn = V, T22 = 0. (b) Lower 
pannels for four MBS: lower left conductance for weak ti2 = 
O.ir (dashed line) and strong ti2 = T (full line); all arc 
equal and only the 1st MBS is coupled to the lead. Lower 
right has the corresponding shot noise. The conductance G 
is presented in units e 2 /h and the shot noise power S in units 
of 2Fe 2 /h. 




FIG. 2: Fano factor F for a double MBS corresponding to 
the upper right panel of Fig.l. The dot-dashed curve shows 
F(V) when only the 1st MBS is coupled to the lead, while 
the solid line stands for case when both MBS couple to the 
lead: fn = F and F 22 = T/2. 



For a given voltage the contribution to the noise is from 
pairs with e™ ax > eV, resulting in jumps at e„ = eV. We 
estimate the averaged current and shot noise using our 
formulas ()16|17[) . The relevant frequencies u> < e™ ax are 
selected by a Fermi function / = 1/(1 + cxp[(w — e n )/S\ 
with small S and we approximate G — [ui — eV + i^'] -1 
with a small 6'. We perform calculations for n < 6. 
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FIG. 3: Conductance (dashed curve) and shot noise (solid 
line) for tunneling into a homogeneous infinite Majorana 
chain as function of U = eV/Y for t = 0.1I\ 



responding SN junction. Wo calculate the Fano factor 
and showed that its value crucially depends on whether 
even or odd number of MBS are involved. The Fano fac- 
tor is similar in magnitude to the SN case, in particular 
F(V — > 0) = 2 for even number of MBS and if only one 
MBS is coupled to the lead, as for SN junction; however, 
when more than one MBS is coupled to the lead (Fig. 2) 
F(0) is non-integer. 
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FIG. 4: The shot noise power as function of a double loga- 
rithmic voltage. The behavior of the noise power has similar 
sharp variations as the current [l5| . 



The shot noise power as a function of a double logarith- 
mic voltage dependence is presented by Fig. 3. S n and 
the average current J show similar behavior, though, the 
jumps arc of different heights. 

Conclusion: We apply the standard Keldysh technique 
to calculate the shot noise power of a Majorana chain 
interacting with a normal lead. We find that the shot 
noise has markedly distinct forms for different numbers 
of MBS and for long ordered or disordered chains. We 
show that for weak coupling to the lead the shot noise as 
well as the conductance are much stronger then the cor- 
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